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Context. The attachment of free electrons to polycondensed aromatic ring molecules (PAHs) is studied for a variety 
of such molecules with different numbers of condensed rings and over a broad range of electron temperatures, using 
a multichannel quantum approach. The calculations of the relevant cross sections are used in turn to obtain the 
corresponding attachment rates for each of the systems under study, and these rates are parametrized as a function of 
temperature using a commonly employed expression for two-body processes. 

Aims. The scope of this work is to establish from first principles the influence of their chemical properties on the 
efficiency of electron-attachment process for PAHs. 

Methods. Quantum multichannel scattering methods are employed to generate the relevant cross sections and hence 
the attachment rates. 

Results. The rates obtained for the present molecules are found to markedly vary within the test ensemble of the present 
work and to be smaller than the values currently used for the full class of PAHs when models of their evolutions in ISM 
environments are employed. The effects of such differences on the evolutions of chemical networks that include both 
PAH and PAH- species are analysed in some details and related to previous calculations. 
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1. Introduction 

The widespread presence in the interstellar medium (ISM) 
of polycyclic aromatic hydrocarbons (PAHs), organic 
molecules made up of several aromatic rings fused to- 
gether, as well as the existence in the same environ- 
ments of their dehydrogenated, ionized and/or protonated 
counterparts, has been largely inferred from the exten- 
sive presence of unidentified infrared emission bands that 
have been observed in the range of 3 to 14 microme- 
ters (iD'Hendecourt fc Ehrenfreundlll997l: iRhee et al.l[2007t 
iParker et alJl2012t iRicks et al.ll2009Tl . Current astrochemi- 
cal models for the PAH formation have largely been derived 
from the combustion chemistry community that suggested 
how these formation processes should occur in warm and 
dense circumstellar envelopes of carbon-rich stars through 
sequential reactions of smalle r radicals with acetylen e at 
temperatures around 1000 K (jCherchneff et all (|l992tl and 
erratum) . Even more recently, their formation at tempera- 
tures down to 10 K has been surmised after an interesting 
investigation combining cross m olecular beam stud ies and 
quantum chemistry calculations (|Parker et al.ll2012n . 

Another low-energy path to possible formations of PAHs 
has come from investigating the role of the free electrons 
which are present in the circumstellar and interstellar me- 
dia a nd in protoplanetary atmospheres (jHerbst fe Osamural 
2008) in order to assess from calculations and possible ex- 
perimental data the likely attachment efficiency of these 
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free projectiles to the PAHs and the ensuing formation of 
very reactive anionic species, which can in turn react with 
radic als present in those environments (jDemarais et al.l 
[20T1 . That even temporary anions of smaller members 
of this family of aromatic molecules could be involved in 
subsequent reactions with cationic radicals has been sug- 
gested also in our earlier computational studie s on benzyne 
anions (|Carelli et al.ll2010t ICarelli et al.ll201lH . Recent ex- 
perimental work has shown that indeed anionic PAHs like 
phenide (CgH^), naphthalenide (CioH^~), and anthracenidc 
(CuHg) can efficiently react with H, H2 an d D2 as observed 
in flowing afterglow-se lected experiments (jDemarais et al.l 
120121 : lYang et al.ll201lt) . It therefore follows from the above 
that to establish the possible effects on the chemistry of 
dark circumstellar regions of the presence of anionic PAHs 
and to model both the possible efficiency of formations and 
the consequences of such additional partners on the evolu- 
tionary history of that chemistry becomes of direct interest 
for improving our understanding of the PAH chemistry in 
the ISM in general. In particular, the present work wishes 
to analyse whether or not the chemical properties of the 
individual component partners in the large series of pos- 
tulated PAHs play any significant role in the evolutionary 
studies of the chemical networks and if it is still a realistic 
choice to view their efficiency as given by a single rate value 
with a unique temperature dependence. 

In the following we therefore intend to approach the 
problem of identifying, if at all possible, a detectable dif- 
ference of behaviour, for a subset of PAHs chosen as initial 
examples for that entire class of molecules, which are re- 
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lated to their differences in chemical properties. In the next 
section we therefore summarize the method we employ for 
obtaining the quantum scattering cross sections for electron 
collisions with the PAHs and also how we use such scatter- 
ing observables to yield the relevant attachment rates. The 
next section III reports our results for a small group of rel- 
evant molecules and discuss their differences and relative 
efficiencies, while section IV describes the final effects of 
such differences on the evolutionary history of the chemical 
networks currently employed for simulations. Our present 
conclusions are finally given in section V. 

2. Computing cross sections and attachment rates 

2.1. Quantum dynamical treatment: an outline 

Both the target one-electron bound orbitals and the im- 
pinging electron wavefunction are expanded as linear com- 
bination of symmetry-adapted spherical harmonics 



A''" 



1 



Ih 



(i) 



where the xfh(@i 0) indicate the linear combinations of 
spherical harmonics Yj 1 and the \p/j.) labels record the spe- 
cific Irreducible Representation (IR) indices [13]. The as- 
sumption that the target molecule is adequately described 
by its ground-state Slater determinant (static exchange ap- 
proximation) leads to the following simplified form for the 
scattering equations: 



V 2 + (E - e) 



F(r) = J dr'V{v,v')F(v') 



(2) 



that can be solved (once the model interaction potential 
for our system is assembled) by using the Schwinger varia- 
tion al method as described in ou r earlier publications ,e.g. 
see (jLucchese &: Gianturcdfl~996l) . The static exchange ap- 
proximation does not include the response of the target to 
the impinging electron, i.e. the correlation and polarization 
effects acting at short and at large electron-target distances 
respectively. That additional part of the overall interaction 
has been modelled in our study by writing 



V cp (r) = | 



ycorr 
ypo 



(r) 
r) 



for r < r match 
for r > r match • 



(3) 



In the present work we have used the Lee- Yang-Parr form 
for V corr (r) [15] and further employed the V pol (r) as a 
function of the polarizability tensor to describe success- 
fully, via the V cp model of Eqn.©, this additional part 
of the overall interaction potential. We further replace the 
exact nonlocal exchange interaction of Eqn.([5]) with a local 
energy-dependent potential. 
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where kp is the Fermi momentum and rj the neutral 
molecule ionization potenti al. We have already found 
(jLucchese fc Gianturcd[l996h that Eqn.(j4]) can yield a rea- 
sonably realistic description of such forces in a great va- 
riety of molecular systems. We can rewrite the scatter- 
ing equations in homogeneous form to obtain the static- 
model-exchange-correlation-polarization (SMECP) approx- 
imation for the scattering event 



i V 2 _ l fc 2 + yst + yep + yFEGE 



F(r)=0. (5) 



V"(r) is the sum of the three local potentials V(r) = 
V st (r) + V cp (r) + V FEGE (r), of which V st (r) describes the 
static interaction of the scattered electron with the tar- 
get's nuclei and electrons, V cp (r) contains the correlation- 
polarization contributions of Eqn.Q and the V FEGE (r) is 
the model exchange potential of Eqn.Q. After integrating 
over the angular variables, Eqn.([S]) takes the form 



1(1 + 1) 



dr 2 



i'h' 



(6) 



and the potential coupling elements are given by 

V iZ'h'( r ) = (x^(r)\V(r)\x7 h '(r)) = 
dr X ^{f)V{v) X P Mr)- 



(7) 



The S-matrix elements, related to the K-matrix, yield the 
set of partialwave phase-shifts 



Si = Si(k) = e 



2iS t (k) 



(8) 



for all the contributing angular momenta. For a scatter- 
ing process characterized by many broad and overlapping 
resonances, the usual Breit-Wigner formula is no longer ef- 
ficient. However it is possible to extract the positions and 
the widths of closely spaced resonant states by using the 
Q-matrix formalism 

Q(E) = inS^ = -ih^ = QHE). (9) 

Eigenvalues of the Q matrix are equal t o the delay time of 
the outgoing wavepacket, as defined in (jSmithlll960f ). For a 
multichannel scattering process the trace of the Q matrix 
is further related to the eigenphase sum 8 sum : 



2 fc^gp = TrQ(E) 



(10) 



where the eigenphase sum is defined via the elements of 
Eqii.©, with k 2 = 2E 



^ sum — ^ $l(E) 



(11) 



1=1 



It can be shown that the first derivative of the eigenphase 
sum can be recast as a sum of Lorentzian functions, each 
associated to a pole of the S matrix, characterized by an 
energy position $l(E) and a width Q{E). This is the proce- 
dure we have followed to analyse the resonance regions in 
our present study. In the present analysis, in fact, we limit 
our study to considering the target molecules to remain in 
the equilibrium geometries of their ground electronic states. 
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2.2. Anion formation rates 

The next step in the calculations involves the thermal aver- 
age over the internal states of the species in the two-body 
(2B) reactions that will be analysed in a moment. By as- 
suming that the free electrons have a kinetic temperature T e 
and that the molecules interacting with them have an inter- 
nal distribution over states which can be associated with a 
temperatur43?M, the existence of local thermal equilibrium 
(LTE) implies, given the large differences between the elec- 
tron's mass and those of the molecules, that Tm ^ T e = T 
(i.e. the molecules are considered as still during collisions) 
and that the rate of formation of molecular anions is given 
by 

kf{T) = g^|^E e ~ e/fcBT / va n (v)f(v,T)dv, (12) 

where a n (v), a function of the relative velocity, is the in- 
tegral cross-section (ICS) for electron-molecule collision la- 
belled by n, a global index that specifies the electronic, vi- 
brational and rotational initial molecular states and the Q's 
are their corresponding partition functions. In the present 
study we limit the target description to be sufficient to gen- 
erate ICS which are summed over all final rotational states 
accessible at the considered energy (and averaged over all 
initial states as taken to be |jo) = 0). The molecule is also 
kept in its vibrational and electronic ground states given 
the likely conditions in the ISM environment. If the f(v,T) 
is taken to be the standard Maxwell-Boltzmann electron's 
velocity distribution, then we can write 

(13) 

which describes the rate of anionic formation via 2B pro- 
cesses and with molecules in their ground electronic, rotovi- 
brational states. 

The actual molecular processes which are possible via 
binary collisions, could be more specifically described as 
those leading first to anionic metastable species formation 
before energy redistribution, i.e. leading to transient nega- 
tive ion (TNI) formation 

M eq +e" i™i>[M eq "]*, (14) 

which corresponds, at specific energy windows, to a reso- 
nant process of attachment before molecular relaxation into 
a bound anion M~ . The dynamical channels of that collision 
process also involve direct scattering event occurring over 
the whole energy range. The competing processes follow- 
ing Ean. (|T4"[l are given by: (i) radiative stabilization (RS) 
dynamics in the chiefly collisionless ISM environment 

[M eq -]* + e" ^ M e ,r + hv (15) 

and (ii) the autodetachment (AD) process before anionic 
stabilization 

[M cq -] ^ M + e~ . (16) 

Current assumptions consider /cad 3> ^rs (|HerbslJ 

1 1 8 lh . which allows us to simplify the expression for the 
overall formation rate kf(T) of Eqn. p^)) 

= / fcTNI fcRS \ 
J V fc AD+fcRS/ 

1 Not to be confused with the kinetic molecular temperature. 



Table 1. Fit coefficients for the moelcules presented in this 
paper, where a is in cm 3 s _1 , (3 is dimensionless, and 7 in 
K. See text and Eqn. (|2TJ)) for further details. 





a 


P 


7 


C6H4 


10.02 


-0.33 


6.69 


C6H5 


1.44 


-0.12 


-14.77 


C14H10 


1.31 


0.21 


-1.12 


C20H12 


3.06 


-0.01 


2.77 


C24H12 


2.74 


0.11 


-1.12 



into the form 

*' m -(£)*-■ (18) 

The above relation then tells us that the /ctni rates can 
provide an upper bound to the true kf(T) insofar as the 
radiative stabilization remains less efficient than the au- 
todetachment process, so that we can safely assume that 

k f (T)<k Tm (T), (19) 

which indicates that the calculations of oo(T) in Eqn. dTBI 
could be approximated by evaluating the whole range of 
elastic ICS plus contributions to TNIs formation at spe- 
cific energy windows. The latter set of cross sections have 
been computed in the present study and will serve us as 
an estimate of the anionic formation rates for PAHs in the 
circumstellar environments. 

3. The computed attachment rates 

In the present, exploratory, analysis we have considered a 
series of molecules which can help us to establish the extent 
of chemical diversity existing among possible PAH candi- 
dates. In particular, we have chosen two single-ring species: 
benzyne (C6H4) and phenil (CeH 5 ) because of their having 
a permanent dipole moment. Thus, although they are often 
consi dered precursors molecules to the PAHs (W eaver et al.l 
120071 ). we analyse here their rates chiefly to test the rele- 
vance of their being polar molecular targets. 

As examples of polyconjugated systems, we report 
here the attachment rates computed for anthracene 
(three rings: C14H10), perylene (five rings: C20H12), 
and coronene (seven rings: C14H12). The full details of 
the cross section calcul ations have been already pub- 
lished for the coronene (ICarelli et all 1201 lh. and are m 
course of pu blication for anthracene (jSanz et all [2012) 
and perylene (|Carelli fc Gianturcoll2012al) . For the single- 
ring m olecules, the benzy ne data have already been pub- 
lished (jCarelli et al.l [201O ) while the phenil cross-sections 
are be ing prepared for publication (jCarelli fc Gianturcol 
2012Q). All the calculations have followed the computa- 
tional method described briefly in the previous Section [2] 

For the 2B reactions, t he expression for the rate coef- 
ficient is usually - e.g. see IWoodall et al.l (|2007l ) - written 
down in the following form 

kf = a U^k) exp(- 7 /T) cmV 1 (20) 

so that the necessary information for a specific reaction 
is stored in the (a,/3, 7) set of parameters identifying the 
temperature dependence of its rate coefficient. 
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Fig. 1. Computed electron capture rates for phenil and 
benzyne polar aromatic molecules. The corresponding fit- 
ting parameters for Eqn. (|20p are also given in the Figure. 



The calculations of FigQ] report, as specific examples, 
the electron-attachment rates for phenil and benzyne, over 
a significant range of temperatures: the parameters from 
the fitting suggested by Ean. ([20| are also given for each 
system (see TabQ]). 

The calculated dip ole moment of C6H4 is f.38 D 
Kraka fc Cremer|ll993l) . while the one of C 6 H 5 is 0.9 D 
McMahon et al.l 20031 ): all the remaining polycondensed 



examples of the present work have no dipole moment. What 
one clearly sees from the data of FigJT] is the strong in- 
crease of both rates at the threshold temperatures: both 
polar targets indicate more efficient attachment for very 
slow electrons, since both systems have been shown to form 
virtual states at zero energy and also dipole-bound states 
for the extra electron attach ed at very low collision energy 
(jCarelli fc Gianturcoll2012~bl ). Furthermore, we see that the 
rates for such PAH precursors vary markedly among them- 
selves and follow the size of the permanent dipole. The lat- 
ter is larger for benzyne and thus the corresponding rates 
for this molecule are a factor of two larger that those for 
phenil at the energy thresholds. Both rates turn out to be 
of the order of 10 -8 cm 3 s _1 in the range of temperatures 
of relevance for ISM conditions (i.e. up to about 100 K), 
although we see that they exhibit differences in their In- 
dependence parameters that are sufficiently noticeable to 
suggest that to use a single choice of a rate for both of 
them would not be accurate enough. In any event, both 
polar molecules indicate strong increases at vanishing en- 
ergies for the existence of enhanced cross section in that 
range of collision energies due to the presence of scattering 
from a dipole potential. 

The situation changes quite markedly when we now 
move to the examples of the polycyclic hydrocarbons with- 
out permanent dipoles which are presented by the rates of 
FigH 

The following comments could be made from the data 
reported in the Figure: 

(i) all three examples show a different temperature depen- 
dence from those shown in Fig [I] since we now see that 
the rates are smallest at vanishing energies and increase 
with increasing T. 

(ii) the sizes of all cross-sections are around 10~ 9 cm 3 s _1 , 
i.e. much smaller that those for polar molecules even if 
the number of C atoms is now larger; 

(iii) the fitting parameters change along the series but the a 
value, somewhat proportional to the rate size, does not 
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Fig. 2. Computed rates of electron attachment to three 
different polycondensed aromatic ring molecules. Top 
panel: anthracene; middle panel: perylene; bottom panel: 
coroncne. See main text for further details. 



increase with increasing the number of C atoms: all the 
cross sections are fairly similar in size and behaviour 
along the small series of PAHs shown in the Figure. 

The present calculations seem to indicate that the chem- 
ical features of the involved PAHs do bear on the efficiency 
of the electron attachment rates obtained from the mod- 
elling that we have reported in the present work. The pres- 
ence of a permanent dipole moment, in fact, is a significant 
rate-enhancer, especially at the lower temperatures, where 
non-polar PAHs are found to have much smaller rates. 
Furthermore, we see that all our computed rates, although 
considered to be upper bound to the true rates for such sys- 
tems, are consistently smaller than the estimates employed 
in ea rlier work for such quantities (|Wakelam fc Herbstj 
I2008h . 

It is therefore interesting to see what would be the ef- 
fect of the present results on the gas-phase models for the 
chemistry of dense interstellar clouds which can now em- 
ploy attachment rates from ab initio calculations. This is 
the aspect of the present work that we shall discuss in the 
following Section. 



4. Testing the evolutionary models 

To now assess the effects of the new rates for the electron 
attachment processes on different PAHs on the evolution 
of the corresponding abundances in the ISM environments 
we have implemented a pseu do time-dependent dense c loud 
model based on the work of Wakcla m fc Herb st ( 2008]) fol- 
lowing their parameters and their recommendations for all 
the other rates involved in the fuller chemical network. We 
discuss below the comparison between the results for the 
overall chemical evolutions using the single value of the rate 
coefficient that they assume for the PAH electron attach- 
ment process (jWakelam fc Herbstl [2008). and the results 
obtained by employing the three rates proposed in this pa- 
per and individually obtained via ab initio calculations that 
we discussed in the previous sections. 
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For the rate coefficient &ea WH2008 suggest two val- 
ues already found in literature that are kf = 10 -7 (iVc) 3 / 4 
cm 3 /s (10montlll986l) and k fa = 1.2 x 10" 7 (A r c ) 1/2 cm 3 /s 



( All amandola et all 



1989), 



where in both cases Nq is the 
number of carbon atoms that form the PAH. For their 
model they assume Nq — 30 a nd arbitrarily d e cide t o use 
the first formula as indicated in lLe Page et al.l (|2001[) . 

What we are therefore presenting below is a test aimed 
at determining the influence of using more chemically accu- 
rate electron attachment rate coefficients on the evolution 
of the ISM model we have employed. In particular, we track 
the evolution of the ISM (i) first with the aforementioned 
fcf rate coefficient and then with the same set of other rates 
but changing that of the PAHs for (ii) anthracene fcfA, (hi) 
perylene kfp, and (iv) coronene kfc- All the other model 
parameters we adopt are the same as those of WH2008, 
namely nn2 = 10 4 cm" 3 , a gas temperature of T = 10 K, 
a cosmic-ray ionization rate of £cr = 1.3 X 10 -17 , a visual 
extinction A v = 10 s , a PAH fractional abundance of 
3.07 x 10 -7 , and an average PAH size of 4 A. We have se- 
lected to show the calculations for the initial species abun- 
dances given by the model EA2 (see Tab.l in WH2008) 
which is an high-metallicity scenario based on the observa- 
tions of the dense cloud £ Ophiuchi, except for He which 
is set to riHc = 0.09 compared to the total hydrogen abun- 
dance. More details of the full choice of parameters can be 
found in WH2008. 

The results of our simulations are displayed in Fig|3]and 
FigHl The first one describes the evolution of the ionization 
fraction as the sum of the number densities of the negative 
species at a given time. WH2008 had shown a fast decrease 
of the total negative charge until it reaches a steady-state 
after 10 3 yr, while the present models have a slower decrease 
in time and reach the steady-state over a much longer time 
scale, i.e. in 10 6 yr or more. The effect of the new calcu- 
lations is also shown in the figure and is rather clear: in 
the present model the PAHs manage to "soak" and cap- 
ture much fewer electrons whe n compared to the results 
from lWakelam fe Herbstl (|2008h . since their selected single 
value for the attachment rate coefficient is much larger and 
thus indicate a more effective process. The PAH - formed 
in this process rapidly recombine with the ions, thereby ob- 
taining a more neutral gas as shown in FigO On the other 
hand, when we replace the WH2008 rate coefficient with 
our separate rates for different molecules we find that the 
neutralization effect is less important, thus obtaining a less 
neutral gas once the steady-state is reached. 

In a similar fashion, the data in FigUshow the evolution 
of the PAH abundances over the time range t = [10 2 , 10 5 ] yr 
and, even in this case, the effect of the new rate coefficients 
is evident. We obtain a larger separation between the abun- 
dances of PAH and PAH - , a quantity that obviously in- 
creases when the attachment rate coeffi cients decrease with 
respec t to the single value chosen in IWakelam fc Herbstl 
(2008]) . 



5. Conclusions 

In the present work we have endeavoured to analyse in 
some detail how one can realistically generate electron- 
attachment cross sections over a broad range of energies 
for gas-phase polycondensed aromatic molecules, a specific 
subset of such molecules chosen to be representative of the 
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Fig. 3. The evolution of the total abundances of the nega- 
tive charged species with time normalized to the total hy- 
drogen number density (i.e. 2 x 10 4 cm -3 ). We compare the 
rates from the original WH2008 model that use their selec- 
tion of kf electron attachment rates (solid) with the cases in 
which use instead: the anthracene rate (long-dashed), the 
perylene rate (dashed), and the coronene rate (dotted). See 
text for further details. 



more general behaviour of PAHs in the interstellar and cir- 
cumstellar environments. In particular, we have selected 
a dynamical model that employs the above cross sections 
to obtain attachment rates over a range of temperatures 
which are of direct interest in the ISM environments. The 
corresponding rates are essentially considered to be upper 
bounds to the true attachment rates, but are nevertheless 
found to be already larger th an those employed in curre nt 
evolutionary models - e.g. see Wakclam fc Herbstl (|2008h . 

Furthermore, the parametric representation of com- 
puted rates is chosen to be the one usually employed 
i n the literature to describe 2B processes in the ISM 
(jWoodall et al.|[2007h and is found by our present work to 
be markedly dependent on at least two chemical properties 
of the aromatic molecules: (i) the presence of a permanent 
dipole moment, a features which increases the size of the 
attachment rates, especially at the threshold temperatures, 
and (ii) only very slowly on the number of C atoms or on 
the number of condensed rings in the molecular system. 

The captured rates have been further employed to 
model the time evolution within a pseudo time-dependent 
de nse cloud descriptio n that follows the one given earlier 
bv IWakelam fc Herbstl (|2008tl . In particular, we have tried 
to analyse the effects of the attachment rate values pro- 
duced by the present work on the evolutions of negative 
charged species with time (e.g. see Fig|3|) and the evolu- 
tions of PAH/PAH - abundances along the same model, as 
well as the electron number density changes over the same 
time period (see FigE]). 

The corresponding calculations suggest rather clearly 
that: 

(i) the attachment rates are smaller, at least over the rele- 
vant range of temperatures, than the values previously 
employed; 
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Fig. 4. Left panel: time evolution of the abundances of PAH (grey) and PAH - (black) normalized to the total hydrogen 
number density (i.e. 2 x 10 4 cm -3 ). We compare the WH2008 electron attachment rates (solid) with those obtained 
using instead: the anthracene rates (long-dashed), the perylene rates (dashed), and the coronene rate (dotted). Note that 
the lines that represent the PAH densities overlap with one another, except those from WH2008. Right panel: electron 
number density time evolution as given by the same set of models. Note that the two panels have different vertical scales. 



(ii) they depend on the chemical properties of the specific 
PAH being considered, especially when comparing an- 
thracene and coronene; 

(iii) the electron densities remain larger over a much broader 
range of time a nd reach a steady state m uch later that 
thought before (jWakelam fc Herbstj[2008l) : 

(iv) the "soaking-up" power of PAHs we have considered 
here turns out to be smaller than expected, hence re- 
ducing the efficiency of the neutralization processes. 

Although we feel that further numerical experiments are 
needed in order to also innclude directly the efficiency of 
polar PAHs within the evolutionary model, the set of ex- 
emplary molecules presented here are already making clear 
to us that chemical variety plays a role which is bigger than 
previously expected and that a molecule-specific modelling 
of the role of electrons in evolutionary studies should be 
sought as much as possible. 
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